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Про визначення початкового наближення 
для контактної задачі 


(Представлено академіком НАН України В. І. Старостенком) 


Запропоновано аналітичні моделі гравітаційного поля і горизонтально-шаруватого гео- 
логічного середовища. Перша модель отримана з плоского нормального потенціалу сили 
тяжіння в локальній точці. Інша модель -- клас Чорного контактних повердонь -- 
отримана з відомого класу Страхова. Подано нові числові алгоритми для визначення 
початкового наблигісення густинного контакту та його асимтптот, у цьому класі. 


Сучасні методи обробки геофізичних даних потребують створення цифрових аналітичних 
моделей поля і геологічного середовища, які націлені на застосування в ГІС. Для цього 
слід виробити відповідні аналітичні конструкції. Одну з них -- модель середовища типу 
контакту пропонуємо до розгляду. 

Визначення контакту у відомих моделях середовища Страхова, Нумерова, Маловичка Ц| 
означає обчислення різниці ундуляцій теоретичної моделі середовища та деякого опорного 
її елемента, від якого ці ундуляції не дуже ухиляються. Але при цьому слід заздалегідь 
знати елементи геометрії (глибину, форму, орієнтацію) опорного елемента (моделі тіла) 
як початкового наближення ітерацій. Для горизонтально-шаруватого середовища виведено 
аналітичні моделі поля і геологічного середовища, в яких ці опорні елементи розраховують 
у процесі ітерацій. 

Аналітична модель поля і середовища. Розглядаємо горизонтально-шаруватий 
простір Е (3) з двома густинними межами дС" й дС» розділу шарів С й С», відстань між 
якими 4і85(0С1, С») 2 0 при 201 (| С» - С». В області Сб - (61,62, 83) - точки облас- 
ті С з С1|)С», С» - С»| )9С», х - (21, 12923) - точки замкнутої необмеженої області 
СЕС, дО о да аві | | ооо 

5 1 дстань між ними дорівнює 7(, Є) ( У (аа - бі) ) ; 
і-і 

Якщо область С» заповнена масами з густиною до (Є), а область Со - Сі «С» -- з гус- 
тиною 64 (Є) (рис. 1), де п(х) -- одинична нормаль до поверхні ДС", зовнішньої щодо С", 
то нормальна, сила тяжіння дорівнює 


д9У(х) д 1 
бта) 1 Ї Юта) б 











4Є, (1) 


3) є 


де 6(2) з 6г(5) - 61(5), Є є С»; (Є) - 0, Є є бо; 7 - гравітаційна стала. 


Силу тяжіння описує функція 


ща) « Інад Уа) « / іі 5(Є вгад, у" Нам Є 4Є, (2) 
4 


оскільки значення сили тяжіння -- це значення модуля градієнта потенціалу |З|, а не її 
складові, в тому числі й (1). Але для малої області О можна використати вираз (1). 


О Ю.І. Дубовенко, 2013 
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Рис. 1. Модель контактного середовища в класі Чорного 


У цій області 0, з відомою мірою наближення |З| ділянку дО - 0 (Д|дС" межі ОСІ" 
вважаємо "необмеженою" площиною: дС - 12: 23 з 23(21, 2) -- соп55, (21,22) Є 50), де 
Зо «з ОПОСЯ. Нехай в околі 07 функції 61(4) - соп8б, Є Є С» С 07, 62(Є) « сопі, 
є«Є0 ПСі,іб -дбг-діа С» -- власна підмножина 0 . Межу 0 обираємо достатньо 
гладкою, для якої 0 з 19: я з 24534, 96), (25, 2) Є 004, 5,7, - 1,2,3; |) 29: --0С32; зі є 
є С (95))). 


За таких умов стосовно середовища викладених вище, отримуємо потенціал тяжіння: 


дУ(Жі2з) | д асуаєзаєз 
осв й 


2123) 


Нерідко маси, що зосереджені в 0 , генерують аномалії и(ї;), які витягнуті, скажімо, 
уздовж осі Ол». Їх доцільно вимірювати по профілях, паралельних О6;. Формула (3) не- 
придатна для обробки таких вимірів. Але, інтегруючи вираз (3) за Є» при умові, що 23 - 0, 
23 з 23(61), а « Є; Є Ь, отримаємо аналітичну модель поля: 


ь Ше Ь 
263 (р - 1)? реа) 
0 о Рак, Фо іа а Го ПЗ В о 4 
це, 0) Щі бу сооучої 23 Щі ВЕРАСИ ОГА З (4) 


ле єФ и (б) і є) 2. 2 (0), ас єр є-- криві, що обмежують знизу і згори тяжіючу 
плоску область 5. 

Аналітична модель середовища є підкласом одного з класів контактних поверхонь, для 
яких гарантовано єдиність розв'язку обернених задач гравіметрії |4-6Ї. 

У класі Остромогильського Оз») (а, Б) однозначно визначають змінну густину (81, 63) 
і контакт 23(61) за значеннями зовнішнього поля |4| для заданого контакту: 


-В-да ; а, 5, 
забрдоа -баз(Є), Є Є (а, 0) (5) 
-й, Сі Я Га, Б, 
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де б23(а) - да3(Ь) - 0, 23(621) -- кусково-неперервна функція; 
|баз(а1)| 2 (ат - а)",  а«атбачо,  а»б, 


|да (21)| 2 с(Ь Пи зано ь-д «алі к« Ь, ан). 


З підкласів осн й осн 


У класі Страхова 5619 (а, Б) |5| контакт задається рівнянням 


легко утворити моделі антикліналей і синкліналей. 


-Тло В ба3(861), 1 Є (а, б), 
23(61) з й ам | - (7) 
-12, 17-90, 


де даз(а) - 0, Ві, По » 0, 23(61) « 0 -- кусково-неперервна функція, задовольняє на кінцях 
інтервалу умовам 
Ібла (Єт)| 2 (Єр - а) "9, | а«фбкачо, оц а», і 
5 
аз(а)| 2-6), б-баакь б», 


набуваючи на (а, б| скінченне число екстремумів. Контакт можна обмежити так: 


за сер ср) Вир вч хе), -25 Пи, зол) Зо с АВ, -2і «Ло, (9) 


а, а, 


де сі, с» Є (а, Р, а точки сі, с» - внутрішні точки відрізку |а, |, або збігаються з його 
межами. У останньому випадку за умови Й « По маємо сі - а, с» - Б. 


Якщо а - -о00, р - Зоо, вирази (7), (8) набувають вигляду 
що 2) і -Т в 622), З є (--Єо, бо); 
|2(Є) -- Ти| Є аз|ЄЇ/?, Є 2 бо, (10) 


|2(6) - Ао| « ва|Є| 8, В.Й; 


де Ві,Пг » 0, а 25 з 2(со) - 8ирі ехіт 2(61)), 2; з 2(са) з ів ехіг 2(61)) при сі, с» Є 
Ї -бо»б0) 1-60 60) 


є | -Єо, 6о|, і (для певності) нехай -25 Є Пі, -2; Є В», 0 «0 Ві «0 По. Для такого контакту 
шарів сталої густини встановлена єдиність |5. 

Нехай С" -- горизонтальна, смуга, що обмежена згори і знизу прямими 2) - Пі, 2 - Па, 
0 «А «А» «2 Зоо, розділена на два шари зі щільністю 6; й до нескінченно гладкою кривою 
дС", з асимптотою 2 - В, п; « В « По», зірчастою щодо нескінченно віддаленої точки, де 


(8) -М «ке, ш а»о, Єр» бо, «Є)єдС. п1) 


Якщо значення Й і стрибок густини Дд - д» - 0; відомі, то контакт ДС шарів за зовнішнім 
полем відновлюється однозначно і є функцією класу (11), який, в свою чергу, є підкласом 
класу (10). Назвемо клас (11) СЬФР( - со, со) класом Чорного (вітчизняного науковця-гео- 
фізика А.В. Чорного, який заклав основи дослідження функціоналів для контактних по- 
верхонь). Клас (10) названо класом Страхова 58078) - со, со). Очевидно, що правдиві від- 
ношення: 


Озі(9) (а, Б) с 5609(а, 6), || СЬО) (оо, со) С 85 со, оо) 5 8509 (а, 6). | (12) 
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Ця модель відображає складний розподіл мас -- п рудних тіл або контактів, які по- 
парно не перетинаються (див. рис. 1). Головною проблемою є єдиність розв'язку задачі (4) 
в класі (11). Ця єдиність гарантована у вузьких підкласах класів (12). Їх можна окресли- 
ти за наявності потужної апріорної інформації про середовище (11). У разі її відсутності 
лишається сподіватися на відомі наслідки з теорем розділення полів, які є справедливими 
і для моделі (11). 

Алгоритми обчислення контакту. Щоб розв'язати нелінійне інтегральне рівнян- 
ня (4), задавши сталі асимптоти, які належать одній з множин (12) допустимих функцій, 
після ряду перетворень отримаємо таке рівняння: 


(ер ші 2-4 Н? 
-е зх (ер - а1)2 в "секр евро ро В ом 


Розглядаючи його як функціонал 4(Х, 2) від шуканої функції 23(Є21)у лінійному просторі, 
знайдемо приріст ДА(7, 2) - (2,21) - (2, 2), де В(а) - А(Б) - 


6 
22 2 
дно з З проти нО | Пе оно" 


362(2) 462(Е) | ) 2 
у о оо |Й АРТЄ94; 14 
пен пе НОРОК я 
та 02(5) «- 2(6) 4 0(2)Р(2). Після ряду аналітичних перетворень на його основі отримуємо 
кілька способів відновлення контакту 23(61), а « Є « Б за заданим полем 4(21), с с х « 4, 
(а, 6) с |с, | за умови плез|а, Б) «/ певі(с, 4. 
Зокрема, при заданих а, б, Н знаходимо Д.Є) с 2пці(Є) - Хп (2), де 


Ь 
/ даній яра 


2 
З 


зауна Є) 
оо НЕ ( УСІ сон (с4є (15) 
оч "00 файні 


При 2па1(2) 2 2п(2) У й (Є), поз 0,1,2,..., Во(Є2) - 0, та відомому значенні 20(5) : 21(2) - 
в р 


Ь 
ол ЗНО т с РР 
З є)2 -- 220 Нтал (6) Є Із рони з (2), (16) 


Ь 
ен" за(Є) Є 





/ ун" 
сані ЗЕМНЕ СОМ 
й / арена / ее Щі 
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Подібним чином, крім контакту 2 - 2(Є), а « Є « Б, знайдемо "асимптоту" Н - соп5і: 
зіприкивдь заднюни Но 1А0(6), нважанічи, що ак з б бла Но НААН, 
п - 0,1,2,..., 1 Д2о(Є) - 0, ДНо - 0, а « Є « 6, отримуємо аналог виразу (15): 


ь 
2242) 
/ "срарзоварото виь ьи Ге пероні 
ь 
ПА БЕ Са с 
й сотО б 


Існують відповідні аналоги для виразів (16) й (17), що отримані за додаткової умови 


Жпа1(6) 72 о (6) 4 Пи (6), Нічна З Но В Ти» 


(19) 
По(2) - 0, то - 0, Н; з Но, РА ОЇ. 


Вибір початкових наближень. У виразах (15) (18) істотною є проблема вибору по- 
чаткового наближення контакту. Розв'яжемо її таким чином. 
Нехай крива (7), с « х « а має на |с, | кілька екстремумів у точках с; Є (с, И|, і - 
щ- 1,2,..., п. Оскільки б » 0, то кожному екстремуму відповідає однойменний екстремум 
д с «(Є), а « Є « фЬ. Визначимо у множині чеку найбільший максимум 25 Є 4(сі) З 
С, 


б воріехі вок і найменший мінімум 2; З 9(с)) З змен (ам Нехай 25 с (сі) з 
с са 

ач Б ЗОЇ зн ее бонено 2(Є)) - відповідні екстремуми 2 (тут абсциси 93 Й 25 
С, с, 


(ді й 2:) збігаються, що не в усіх випадках вірно. 
Початкове наближення контакту задаємо двома способами: 


хо з 20(2) - -Но КАЄ, Но), ак Є (20) 
20 - 20(2) - -Но ян Ко ДА(2), ак З «є Ь, (21) 
де Но » 0 -- задана стала "асимптота" 20 на лівому кінці (а, |, Д(Є) - (2) - 4(а), 
44(5, Но) - аналітично продовжене значення ДА(Є) на рівні - Но, ко - К(Но) - невідо- 


мий коефіцієнт. 
Знайдемо Ко. Нехай у точці с Є (а, Б| розташована 23 (2;). Отримаємо такі нерівності для 
визначення початкових 25 й 2; при заданих Н, 93 Й 9: 


(1 - су? Н? сон? Н 
20 - сі) і зано о ЕН агсбв зон 
48 2 | і ветаевуюаго зд ЧНО на оо) 
25 
-2 зво Зб --о 22 
СВ р ( ) 
фе ор ав. Н 
«(1-сі)1 сії -2Н агсі - 
з«(170і) (реадбней ба; ні ос) 
я 
-2 ПТ р 23 
г ахсбє рожа разу (23) 
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Відшукання значень 23 й 2; однотипне, тому індекси опущено. Утворимо функцію 


пауз ач Фа з ен неї? а) два оо 


-а-сш - 22 - НУ - сс -- Н2) - ЗНагсів -9. (24) 


Не 1-0) 
Наближення 2, якогось із значень 235 й 2; визначаємо за схемою Чебишева |?|: 


Ж За Ган) - (вад (аа) 
тек т Рв) 2 РІЗ (жи) 3 


де /(2) визначається з функції (24). Процес (25) при вдалому виборі 2, швидко збігається: 
4 ітерації для отримання 2 9) з точністю до 10 7, що узгоджується з результатами праці |7. 

Зауваження 1. З виразів (22), (23) випливає, що /(23) 2 0, (2) Є 0, аза не дуже малих 
Н» 0 маємо /(25).» 0, (2) « 0. Тому корені 2в, єн рівнянь /(2в, 95) - 0, /(2н, ді) - 0 
дають наближення 23 й 2; знизу та зверху, що погіршує "якість" початкового наближення 
2о з- 20(Є). 

Зауваження 2. Величини 23 (2і) визначаються неякісно при великих Н, тим паче, якщо 
вони є великими. Тому разом із 25, 2; визначають при бд - 9(Р) - (а) 7 0, крім заданої 
лівої Н" і праву "асимптоту" Н" кривої 2 -- 2(Є) за описаною вище схемою, тільки у виразі 
10) замість д5 (ді) беруть ду для визначення Н". 


(25) 





Обчисливши для кожного заданого значення Н., 2 - 1,2,...,т, лівої "асимптоти" шу- 
каної кривої 2 - 2(Є), а « Є « Б, значення параметрів 29, 29, НЕ 2--1,2,...,т, знайдемо 
ці) -20-НІ, | щу тан, | АНу-НІ- НУ, 
9 -- 9(6) - (а), || дфтах - 95 - (а), || ддшію 7- 9 - 9(), (26) 
(б в. ФГ су 0 Їщо| АНІ 
звідки о ро) де Кі Ко же Кз за точо. 
«їх і З боупах| на ддааіо | І94| 


Виведено нові аналітичні конструкції для обчислення складних контактів (20), (21). Для 
їх обчислення запропоновано процес простих ітерацій (25), а для послідовного уточнення 
виконують процес (26). 

Схему відшукування початкового наближення для контактів з класів "антикліналей" 


і 12 
і "синкліналей" Ов ща, Ь) апробовано на функціях, які моделюють антикліналь і синклі- 


наль: Є) - -Наа(1 - Є)Є, 220) - -Н - а( - 2)2, 0 « Є « 1. Ці конструкції краще 
збігаються при обчисленні за першим зі способів (15)- (17) порівняно з відомими |8|. Для 
асимптот Н;, Н, Рі характер збіжності ітерацій є аналогічним. 
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Ю.И. Дубовенко 


Об определений начального приближения для контактной задачи 


Предложетнь, аналитическиє модели гравитационного поля ш горизонтально-слоистой, ге0- 
логическої, средь. Первая модель получена из плоского нормального потенциала силь; тя- 
ожсеєтли в локальной точке. Другая модель -- клас Черного контактньх повераностей -- 
получена из известного класса Страхова. Дань, новье численнье алгоритм для определе- 
ния начального приблигісения плотностного контакта и его асимптот, в атом классе. 


Уп. І. рибоуепКо 


Оп бе іпібіа! арргохіпабіоп даебпібіоп бог а сопбасі ргобіена 


Апаіубіса! тодеіз Їотг іпе угауйу Дей апа їВе Рогігопіайу Іауетеа деоіодіса! тейішт, ате обетеа. 
ТРе Ійг5і тодеї і8 детіуеа гот, їйе 9-Р поттаї дтауйу роїспіза! а? а Їоса| роті. ТПе іазі тодеі -- 
а СВотпу сопіасі зит/асез сіавз -- 5 йетічей Їтот, їе Кпоцт ЗігаКНог сіавз. Меш питетіса! адо- 
тіййтз от іПе дептийоп ої пе зіаті арртохітабіотп ої а аепзібу іпіет(асе апа її5 азутріоїез іт, 
іа сіазз ате дічет. 
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